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(57) Abstract 

A special mobile radio mobile, vehicle 
or personal location system is described where 
the system operates to locate any one of the 
transmitters in an SMR system. The present 
invention operates by identifying the received 
signals, removing the modulation information 
to reconstruct a demodulated carrier and ana- 
lyzing the carrier signals for muitipath distor- 
tion. The muitioath disruptions are removed to 
calculate the true incident angle of the trans- 
mitted signal at a base station. A plurality of 
base stations are used to simultaneously cal- 
culate the incident angle of the transmitted 
SMR signal and a central location calculates 
the most probable location of the transmitter 
based upon cri angulation. The incident angle 
of the transmitted SMR signal at each base 
station is performed by calculating the phase 
differences of the incoming voice signal with 
phase ambiguity minimized by using a plural- 
ity of dipole antennas spaced at an irregular 
distance to eliminate phase ambiguity error. 
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DIRECTION FINDING AND MOBIT F LOCATION SYSTEM 
FOR TRUNKED MOBILE RADIO SYSTEMS 



FIFTH OF T HE INVENTION 
The present patent application relates to radio transmitter 
location in a communication system and in particular to locating a transmitter 
10 in a special mobile radio system. 

BACKGROUND OF THE INVENTION 
Mobile trunkiog radio is well known in the art and is also 
sometimes known as Special Mobile Radio (SMR) systems. A trunked radio 
15 system is a system setup to serve a plurality of users in a licensed radio 

frequency band to provide radio or radiotelephone services to base or mobile 
units within a defined geographic market area Typically SMR systems are 
used in the public safety industry for such services as police, fire, and medical 
providers and also in the industrial or land transportation industries, such as, 
20 taxicab services, bus services and the like. Such an SMR system typically 
also includes repeater stations in which the signals from individual users are 
amplified and repeated at various places within the geographic market to 
provide a wider geographic accessibility. 

In the United States, an SMR service provider typically obtains 
25 ^a license for a band of frequencies within the 800 and 900 MHz frequency 
range from Federal Communication System (FCC). The SMR system is 
described as a trunked system which means that several channels are allocated 
within a given frequency band and the transmitters and receivers use 
automatic switching so that all users can access any one of the allocated 
30 channels which are not currendy in use. This form of frequency sharing - 
results in a minimum wait or blockage of a message through the system by an 
individual user since the user's transmitter shares all available channels. 

By way of example, the FCC has allocated 250 channels 
located in the frequency band of 809 to 816 MHz. Each channel spacing is 
35 25 KHz and the maximum deviation of a transmit frequency is plus or minus 
5 KHz. The required frequency stability is 2.5 ppm for transmitters. The 
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FCC allocation also provides for repeater channels which are located 45 MHz 
above these frequencies. For purposes of the present discussion, , repeater 
frequencies are not required and hence will not be discussed 

Typically the mobile transceivers in an SMR system can be 
5 linked to a public switched telephone network by base station or to cellular 
telephone systems. Access to the public switched network may be from a 
variety of means including DTMF tone signalling for automatic operation 
without human intervention. All the SMR channels are voice channels for 
carrying analog voice signals to in-band tone signalling may be used to 
10 interface to a land-based public switched telephone network. Links to cellular 
telephone systems has been done in systems described as Extended Special 
Mobile Radio (ESMR) systems. 

Special mobile radio systems to which the present invention 
applies allow the automatic sharing of a plurality of channels by all users of 
15 the system in a blocking communication model. . The system is blocking since 
there are typically more users than there are channels and the channels are 
allocated based on demand Statistically the number of channels allocated to 
a user group is based upon a peak demand usage to minimize blockage of the 
; channels due to unavailability of an open channel. Of course the greater 
20 number of channels in a special mobile radio system, , the fewer times the 
channels are blocked 

Another feature of the special mobile radio system is the use of 
a control channel to allow coherent channel allocation to maximize usage of 
the system. The control channel may be a dedicated frequency, a subcamer 
25 frequency or a sub-audible frequency in selected voice channels in the SMR 
system Allocation and usage of the communication channels is therefore 
negotiated through the control channel in a distributive processing topology in 
which a central control site and the various user sites negotiate, the use of the 
channels during handshake transmissions prior to channel usage. Various 
30 types of communications protocol may be employed to ensure non- 
interference among the various users of the system. The control channel 
typically carries control information in a digital format. 
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Prior art Figure 1 shows a typical simplified SMR radio system 
used to communicated with a plurality of cars 101a, 101b, 101c, etc. from a 
control station 103 via communications antenna 102. In this system, all the 
mobile radio transceivers 101a, 101b, 101c share a number of common 
5 frequency channels in which an individual transceiver may use any one of the 
unused channels at any given time to communicate with the control station. 

There is a need in the art, however, to locate the mobile 
transceivers in a special mobile radio system. There is a need to identify and 
locate the mobile transmitters in an SMR system while the mobile transmitters 
10 are transmitting. There is also a need in the art to track the movement of the 
transmitters even if the transmitters change frequency channels over a period 
of time. 

Location of radio transmitters is a technology that has been 
used for many years and in many applications. For radio location, die use of 

15 mobile transmitters and a plurality of fixed base station receivers is generally 
known in the art where the fixed receivers locate the source of the radio 
signal by triangulation of the radio signals. Triangulation is a technique 
where the fixed base station receivers each determine the approximate 
direction of the radio signal from the transmitter. The approximate direction 

20 lines are then drawn on a map where the lines will cross. Since the direction 
lines are approximate,, and typically three fixed base station receivers are used, 
the intersecting area between the three lines is usually shaped like a triangle 
since the lines rarely meet at a single point The transmitter is then located 
. within the triangle. 

25 There is a need in the art for a ground-based vehicle or mobile 

transmitter location system which is capable of locating SMR transmitters. 
There is a. further- need for a ground-based mobile transmitter location system 
which is immune from disruption of the radio signals in an urban or irregular 
terrain environment. For these reasons, and for other reasons stated below 

30 which will become apparent to those skilled in the art upon reading and 

understanding the present specification, there is a need in the art for a special 
mobile radio vehicle or personal locator system. 
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SUMMARY OF THE INVENTION 
The present invention solves the above-mentioned problems 
with locating transmitters in special mobile radio systems and other problems 
which will be understood by those skilled in the art upon reading and 
5 understanding the present specification. The present invention describes a 
mobile, personal or vehicle location system which locates any one of a large 
number of transmitters in an SMR system. The present invention operates by 
first identifying the transmitter during the SMR transmissions. If the location 
of the transmitter requires periodic recalculation, the present invention tracks 
10 any channel frequency changes based on that particular transmitter's use of the 
frequency channels. The received signals are analyzed to calculate the 
incident angle of the transmitted signal at a base station. A plurality of base 
stations are used to simultaneously calculate the incident angle of the 
transmitted SMR signal and a central station calculates the most probable 
15 location of the transmitter based upon triangulation 

To perform the triangulation at the central station, the incident 
angle of the transmitted SMR signal is calculated at each base station. The 
received SMR signal is detected, the transmitter identity is determined and the 
modulation information is removed. Once the base carrier signal is recovered, 
20 the phase of the received carrier is compared between a plurality of antennas 
in an irregular spaced array of antennas. The phase differences of the 
incoming carrier signal are compared between the antennas and the phase 
ambiguity is minimized to accurately determine the incident angle of the 
received carrier. 
25 , 

BRIEF DESCRIPTION OF THE DRAWINGS 
In the drawings, where like numerals refer to like components 
throughout the several views, 

Figure 1 is a prior art description of a typical Special Mobile 
30 Radio (SMR) system. 

Figure 2 is a block diagram of the present invention adapted to 
locate SMR transmitters by triangulation 
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Figure 3 is a diagram of the receiving antennas of the base 
stations including a plurality of dipole antennas and a s ing le reference omm 
antenna 

Figure 4 is a block diagram of a carrier recovery circuit 
5 Figure 5 is a three antenna array for an inlerferometric direction 

finding system- 
Figure 6 is an ambiguity plane plot of the three antenna array 

of Figure 5. 

Figure 7 is the ambiguity plane for the ratios 5:6 of A/2 the 
10 three antenna array of Figure 5. 

Figure 8 is diagram showing triangulation of the location of a 
SMR transmitter between three base stations. 

Figure 9 is a diagram of a gnomonic plot of the measured 
bearings of the direction finding algorithm on a sphere. 

15 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
In the following detailed description of the preferred 
embodiment, reference is made to the accompanying drawings which form a 
part hereof, and in which is shown by way of illustration specific preferred 
20 embodiments in which the inventions may be practiced These ^bodiments 
are described in sufficient detail to enable those skilled in the art to practice 
- the invention, and it is to be understood that other embodiments may be 
utilized and that structural, logical and electrical changes may be made 
without departing from the spirit arid scope of the present inventions. The 
25 following detailed description is, therefore, not to be taken in a limiting, sense, 
and the scope of the present inventions is defined only by the appended 
claims. 

. The present invention shown diagrammatically in Figure 2 
operates to locate the mobile transceivers in a special mobile radio system. 
30 The transceivers may be located in a car 101, in a building (not shown), 

carried by a person (not shown), carried within containers or packages, or any 
number or variety of mobile carriers. In a single transmission, the SMR 
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transmitters are transmitting on a single frequency channel within the SMR 
channel space. Over an extended period of time, the transmitter may be 
transmit a plurality of discrete messages using a variety of channel 
frequencies which may be different for each messages. The location 
5 algorithms monitor the control channel and, if necessary, track the channel 
frequency changes in the SMR system if a continuous tracking of the 
transmitter is required 

The preferred embodiment of the present invention is designed 
as a separate but parallel operating system used to locate SMR transmitters 
10 without any modification of an existing SMR system. The present invention 
is designed, in effect, to be an "add-on" system to add the capabilities of SMR 
transmitter location without any modifications to the existing SMR system 
whatsoever. Those skilled in the art will readily appreciate, however, that the 
present invention may be made a partial or integral part of an SMR system to 
15 save on duplication of hardware and installation costs. 

Typically, the location of the transmitter is determined once but 
for a moving vehicle, periodically updated location information may be 
required. Usually, the SMR transmitters are activated only by the operator. 
In a system which is designed to include the present invention as an integral 
20 part of the overall SMR system, however, the SMR transmitters may be 
polled by a base station or central station to transmit without operator 
intervention so the SMR transmitter location may be determined. This would 
be a "fleet location poll". Also, the present invention may be used as a "panic 
button" or "assistance request" system in which SMR transmitters are used 
25 summon help on operator demand, or automatically summon help such as in 
the case of a crash-activated deceleration switch similar to the deceleration 
switches used, to activate vehicle airbag systems. 

As shown in Figures 2, a plurality of base station antennas 
202a, 202b and 202c receive the SMR radio transmissions and the base station 
30 equipment 203a, 203b and 203c, respectively, each locate the incident angle 
of arrival of the true signal even in the face of selective fading and multipara 
receptions. Several of the most probable incident angles of the received 
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signals are sent from the base stations 203a, 203b and 203c to a central site 
204 for triansulatioiL See also Figure 8. 

The location system comprises a plurality of base stations (at 
least two) which transmit signals to the SMR remote mobile units and receive 
5 signals from the SMR remote mobile units. The base stations communicate 
using standard SMR channels carrying voice band information to the SMR 
remote mobile units, or the base stations may simply send alert or 
interrogation commands to activate the SMR remote mobile units. Thus the 
SMR remote mobile units may be activated locally (for example by the 

10 wearer or vehicle) or remotely by the base station. The SMR remote mobile 
units may be used by a person as a communications radio, a radio-telephone, 
a cellular telephone link or as an emergency locator. The SMR remote 
mobile units may also be used to locate the vehicle if stolen, in an accident, 
for vehicle tracking or as an integral communication device. 

15 The SMR remote mobile units communicate over licensed 

channels in an SMR system using a standard SMR infrastructure. SMR 
systems are manufactured by vendors such as Motorola and E. F. Johnson. In 
the preferred embodiment of the present invention, the identification and 
location determination components of the present invention are designed to be 

20 installed as a parallel system operating completely independent of the existing 
SMR system. In this fashion, the present invention may be installed and 
operated independent of the SMR system or it may be made an* integral part 
of an existing SMR system 

The transmitted signals from the SMR remote mobile units are 

25 received by the base stations which are then used to locate the SMR remote 
mobile units. An array of irregularly-spaced dipole receiving antennas is 
added to the existing omni receiving antenna at each base station. For 
direction finding, the existing omni antenna is used as a reference antenna 
This is an optional feature since, as will be described more fully below, any 

30 one of the dipole antennas of the array may be used as a reference antenna. . 
The direction of the received signals from the SMR remote mobile units is 
determined by interferometric direction finding using phase difference 
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comparisons between the various dipoles of the array. Unique algorithms 
using histogram analysis and spatial filtering are performed, if necessary, on 
the received signals at the base stations to eliminate the ambiguity caused by 
signal interference, signal drop out, ambiguity in bearing measurements 

5 received from behind the receiving antennas and multi-path reflections. 

The present invention is used in a wide variety of applications 
wherever a special mobile radio system infrastructure already exists. The 
present invention may be installed on the radio towers already used to serve 
the SMR system without affecting the operation of such a system. The 

10 present invention may be used to locate vehicles equipped with SMR devices, 
persons carrying such radios of SMR radios located in homes or buildings. 
The system may also be used as a "panic button 1 locator where an operator is 
in need of assistance but may not be able to communicate his identity or 
location. 

15 System Overview 

The present invention allows the accurate location of the SMR 
transmitters even in noisy and multi-path environments such as those found in 
urban environments and irregular terrain environments. 

20 Base Station Design 

The analysis and decoding of the signals received by the base 
stations from the SMR remote mobile units is done almost entirely in the 
digital domain The carrier frequencies of the SMR frequencies are down- 
converted to a lower frequency in each base station and are then digitally 

25 sampled. The digital samples are processed to locate the phase information of 
interest for direction finding. Each base station is constructed with a plurality 
of digital signal processors which enable direction finding of a plurality of 
simultaneously transmitting SMR remote mobile units transmitting on 
different frequency channels. 

30 The receiving antennas of the base stations are constructed 

using a plurality of dipole antennas 401a, 401b, though 40 In (generally 
- referred to as 401) arranged in an array, as shown in Figure 3 and a single 
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reference omni antenna 404. The orani antenna 404 is not necessary for the 
preferred embodiment of the present invention since any one of the dipole 
antennas of Figure 3 may act as a reference antenna for any of the other 
antennas. The omni reference antenna 404 is used as a reference antenna 
5 because the omni antenna 404 exists as part of the SMR receiving and 

transmitter infrastructure upon which the present invention relies and is used 
for receiving the SMR voice channel messages from the SMR remote mobile 
units. 

In the preferred embodiment of the present invention, the 

10 direction finding (DF) antennas of each base station uses twenty-four sector 
direction finding dipole antennas 401 arranged in four linear arrays as shown 
diagrammatically in Figure 3. Each of the four linear arrays consists of six 
dipole antennas and an omni reference antenna. The phase between the 
antennas is measured indirectly - each instant of time one of the 24 antennas 

15 is selected and its phase is compared to the omni antenna 

The switching of the antennas is done every 5 milliseconds 
during the transmission of a particular SMR transmitter. The antennas 401 
are continuously scanned since there could be several SMR remote mobile 
units transmitting simultaneously from different directions. Differentiating 

20 between the various simultaneously transmitting SMR remote mobile units is 
possible. As described above, the simultaneously received signals from 
different SMR transmitters will be on separate channels and the base stations 
are capable of distinguishing and tracking the SMR transmitters 
simultaneously. This is due to the use of a plurality of digital signal 

25 processors operating in parallel to analyze the incoming signals using Fast 
Fourier Transforms (200 Hz bins). The direction finding process, described 
more fully below- is implemented in digital signal processors (DSPs) located 
within each base station and a general purpose computer located in the central 
site. 

30 SMR System Design Considerations 

The present invention must be implemented to operate 
according to the frequencies, modulation techniques and control protocol of 
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the specific SMR system with which the present invention will operate. For 
example, the present invention will track the mobile SMR transmitters which 
are using only specific frequency channels, specific modulation techniques and 
specific identification techniques. For purposes of describing the present 
5 invention, but not by way of limitation, the present invention will be 

described to operate on an SMR system which transmits in the 800-900 MHz 
band using FM modulation of analog voice information on the earner. The 
identity of the transmitter is in an audio subcamer on the modulated canier. 
As a first step, the identity of the transmitter is easily 

10 determined by either using the existing SMR system or by duplicating the 
identity determination function of the SMR base station in a separately 
operating component. The existing SMR system places the identity of the 
transmitter on an audio subcamer although those skilled in the art will readily 
recognize that the ID information could be part of a message preamble, a sub- 

15 audible carrier or transmitted in a parallel control channel. Regardless of the 
technique used, the identity of the transmitter is automatically transmitting and 
automatically determined at the receiver site. 

As a second step, before the location algorithm can begin the 
processing of the phase information of the carrier to determine the incident 

20 angle of the received carrier, each base station must detect the presence of the 
carrier so the base stations are synchronized in their operation Each channel 
in the SMR system will naturally contain background noise and the direction 
finding algorithms should not be allowed to "free run" on noise. Instead, the 
carrier signal must be detected and tracked to signal the start of the direction 

25 finding algorithms in each of the base stations. The detection of the carrier in 
an SMR system is quite simple. ■ In the case of FM modulation of the carrier, 
a simple FM detector circuit is used to detect the presence of an FM carrier. 
The detection of the carrier must be done in each channel where the target 
transmitter is transmitting in the case where the target transmitter is changing 

30 channels. 
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Phase Determination of the SMR Carrier 
Phase comparisons of the received carrier are calculated as 
between a plurality of antennas in an array to determine the incident angle of 
the incoming carrier signal at each base station Since the SMR carrier signal 
5 is severely modulated, phase comparisons cannot easily be made on the 
modulated earner. Thus the modulation information on the carrier must be 
removed before accurate phase comparisons are made. By removing the 
modulation information, the earner signal is treated essentially as a sine wave 
in the calculations described below. 

10 Figure 4 is a block diagram of a carrier recovery circuit The 

modulation information is removed and the carrier signal is recovered entirely 
in the digital domain. The modulated SMR carrier signal is received at the 
omni antenna and one of the direction finding antennas and the signals are 
downcorrverted to a lower intermediate frequency such as 75KHz. The 

15 downconverted carriers are digitized by an AID converter at an approximately 
iOOKHz sampling rate. A digital high-pass filter and decimation filters 
further limit the carriers and the conjugate of the DF carrier signal is mixed 
with the omni carrier signal to produce a near-sinusoidal reproduction of the 
unmodulated carrier. 

20 The modulated carrier is also a relatively wide band signal . 

which needs to be narrowed The frequency modulated carrier may have a 
frequency deviation of plus or minus 5 KHz around a center frequency. The 
phase comparisons between the antennas may not be able to track phase 
changes accurately enough if the slope of the phase changes in the modulated 

25 carrier is too severe. Thus the unmodulated carrier is narrow banded by the 
high pass filters as shown in Figure 4. The direction finding algorithms 
described below are generally immune to some phase distortion and the 
narrow banding of the carrier aids in the accuracy of the overall process. 

The SMR carrier is also subject to selective fading, drop outs 

30 and multipath distortion especially if the transmitter is in a rapidly moving 
vehicle. Thus the carrier may not always be available and the direction 
finding algorithms may be required to operate continuously to update the 
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location of the transmitter as it moves. The steps of detection, identification, 
demodulation and direction finding are performed periodically at preset 
intervals and can also be restarted in the face of interrupts. 

The determination of the phase of the carrier must first be done 
5 before a comparison of the phase of the carrier can be made between 

antennas. The detection of the phase information of the carrier is performed 
in the preferred embodiment of the present invention by first comparing the 
received phase of the demodulated carrier at the omni reference antenna to the 
received phase of the demodulated carrier at one of the dipole antennas in the 
10 array. A matching filter or autocorrelation function is performed on the two 
received signals until a reasonable match of the phase of the two received 
carrier signals, is determined This provides the carrier inforcoation which is 
then used to synchronize the phase comparisons between the plurality of 
antennas of the direction finding array . 

15 - * 

Antenna Design at the Base Stations 
The DSPs of the base station system compare the exact phase 
amplitudes of each of the dipole antennas of the four arrays to the phase 
amplitude at the reference antenna. In the antenna array shown in Figure 3, a 
20 circle surrounding the array is divided into eight sectors. The actual direction 
computation is done in the sectors according to Table 1. i 

TABLE 1 : Sector Allocation on the DF Antenna Array 

25 30° to 60° - top and left arrays 

60° to 120° - left array 

120° to 150° - left and bottom array 
: 150° to 210° - bottom array - . 

210° to 240° - bottom and right array 
30 240° to 300° - right array ' 

300° to 330° - right and top array 

330° to 30° - top array 

35 Those skilled in the art will readily recognize that the whole 

direction finding array does not have to be installed in a single locatioa For 
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instance, each linear array of six dipole antennas could be installed on a 
different face of roof of a high rise building. In some cases an SMR remote 
mobile unit cannot transmit from specific sector. For instance one of the 
sectors faces the sea In that case the direction finding array could be 
5 reduced, and only a single face of the direction finding antenna array will be 
installed at that sector. Usually the interferometer, does net operate in the case 
where an SMR remote mobile unit is very close to the base station (less than 
1 kilometer). In such a condition, the other base stations will still perform the 
direction finding and the central site will still be able to perform the 
10 triangulation. 

For the reasons exemplified more fully below, the spacing of 
the six dipole antennas in each of the four sides of the array is critical to 
maintaining a low degree of ambiguity in direction finding and multi-path 
rejection. The base ratios of the base lines was found by simulations with 

15 several relative prime numbers and is selected to be 6 : 10 : 4 : 3 : 5 to 

produce a phase gain PG = 56tl The ratio was computed through extensive 
simulations. Although the spacing between the various antenna is critical to 
accurate determination of the incident angle of the carrier, the order of the 
spacing of the antennas in the array is not critical. The order chosen here was 

20 , based primarily upon a physical implementation consideration that the array 
frame be mechanically balanced 

Interferometric Direction Finding 
The direction finding technique used in the preferred 
25 embodiment of the present invention is performed by measuring the electrical 
phase difference between the antenna of the array caused by the angle of the 
incoming wavefront of the SMR frequency of the channel currently being 
used by an SMR transmitter. In ideal conditions where there is no multi-path 
disruptions of the wavefront the incoming wave from the SMR remote mobile 
30 unit is coming in angle of 9 from the perpendicular line to the base line on 
which the antennas are mounted 
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Figure 5 is an ©cample of how the interferometric direction 
finding technique of the preferred embodiment of the present invention is^~ 
practiced Figure 5 shows, by way of example and not by limitation, only 
three antennas of the array. The three antennas are used to determine phase 
5 differences of the incoming signals between pairs of the antenna- As 
described above, the phase difference may also be measured between each 
antenna and a reference antenna such as the omni antenna. 

The interferometer direction finding (DF) technique calculates 
the direction of arrival of the wavefront by measuring the phase of the 

10 incoming wave front in several places in space. This interferometric method 
is a planar interferometric method which assumes that the antenna array and 
the transmitting SMR remote mobile unit are all in the same relative plane. 
Differences in height of the SMR remote mobile unit and the antenna array 
are ignore and cause little error in the direction finding of the present 

15 invention. 

As seen from Figure 5, the three antenna 501, 502 and 503 are 
irregularly spaced such, that distance between antennas 501 and 502 is 
great than distance D 12 between antennas 502 and 503. The distance between 
the antennas is selected to be a multiple of a prime number sequence to 

20 ensure an irregular spacing of all antennas in the array which increases the 
probability that an electrical phase difference will always be measured 
between some of the antenna members of the array. The distance between the 
antennas is important to minimize ambiguity errors in measuring the phase of 
the incoming wavefront. If the antennas are spaced too widely apart, there is 

25 increased ambiguity in determining whether the phase difference of the same 
portion of the voice-modulated signal is being compared or whether the phase 
of different portion of the voice-modulated signals are being mistakenly 
compared Since the licensed frequencies of the SMR channels operate in the 
800 MHz bandL the wavelength X of the incoming signal is approximately 38 

30 centimeters. 
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Again, referring to Figure 5, the phase difference between the 
two antennas 502 and 503 is know as D 12 and the formula for measuring the 
angle of the incoming signal is : 

5 (1) (cp 12 +27ii: i2 ) /2tz =D 12 SIN(e) /X 

where 

0 is the angle of the incoming wavefront measured from the 
10 perpendicular of the antenna anray baseline, 

(p^ is the electrical phase measured between antennas 502 and 
503 

K 12 is an ambiguity factor since the system is capable of 

measuring (j> between ±.7i since the antennas could be 
15 measuring the phase difference of more than one carrier 

cycle \ 

D i2 is the distance between the antennas 502 and 503, and 
X is die wavelength of the incoming signal. 
According to this formula, as long as IS smaller than A/2, the ambiguity 
20 factor K 12 equals 0. Thus the largest antenna spacing in the array is selected ; 
- to be smaller that 15 centimeters. 

In order to estimate the inherent instrumental error of the 
direction finding interferometer, the above equation can be written in a 
different form: 

25 

( 2 ) SIN(Q) = (cp 12 +27tiC 12 ) / {2nD x2 /X) 

Assuming that all the errors in estimating 0 is due to inaccuracies in 
estimating (p, differentiating the above equation produces: 
30 , 

(3) . a(6)=a{cp 12 )/(27c£? 12 COS{6)/X) 
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where a is the standard deviation. ; 
From this equation, inherent limitations in the system are determined The ! 

i 

electrical error is divided by the factor 2nD 12 /A, which is termed the phase j 

s 

gain (PG) which reduces the enror by the ratio between the distance and the j 

r 

5 wave length. Also, the error rapidly grows to infinity as 6 approaches 90°. 
Thus the coverage of the present linear interferometer is limited to ± 60° as 
shown by the division of the anray into sectors described in above in 
conjunction with Table 1. In order to reduce the error, the phase gain must j 
be made as large as possible and each antenna face of the antenna array of the j 
10 present linear interferometer array covers only 1/4 of the space. The antenna 

r 

spacing in the array must have a distance smaller than A/2 in order to resolve j 
ambiguities. j 

Ambiguities in Interferometer Directi on Finding 
15 There are two types of ambiguities in interferometer DF: front 

to back ambiguities (the linear interferometer only covers effectively 1/3 of its 
surrounding area) and multiple possible solutions of the interferometer 
equations, in case where the distance between two antennas is greater than 
A/2. This latter ambiguity is caused by measuring the phase difference 
20 between different portions of the voice-modulated signal of the incoming 
wavefront 

The front to back ambiguity is solved in the present invention 
by using multiple arrays, as shown in Figure 3. In order to cover 360°, the 
present invention uses four sub-arrays of six dipole antennas each, for 

25 covering 90°. The antennas in each array have a vertical beam width of 120°. 
The base station scans all the antennas of the array through the 1:24 antenna 
selector 402 and selects the proper sub-array by comparing the amplitudes 
received from antennas in each array. If two arrays are almost equal in 
amplitude (within 3 dB), then the interferometer process is done on both 

30 arrays and the direction of arrival (DO A) of the signal is computed either 
from both or from the proper array which is selected according to the 
direction of arrival closer to its perpendicular. See Figure 8. 

i 
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The phase ambiguity is solved by using more than two antennas 
in each linear array to compare the incoming voice-modulated carrier signals. 
By way of example, but not by limitation, a three antenna array is as shown 
in Figure 5. The distances between the antennas that are elected as D l2 =r 2*A/2 
5 and D23=3*A/2. Using these ratios, the equation labeled (2) above becomes: 

(4) SIN{6) =<f> 12 /2n+K 12 

10 (5) SIN(d) =923/3* + 2JC 23 /3 

Figure 6 is an ambiguity plane plot of the three antenna array 
of Figure 5. Figure 6 shows the phase difference (fl 2 plotted against the phase 
difference (fo as an interferometer ambiguity plane plot Figure 6 shows a 
15 family of lines, .601, 602, 603, 604, and 605, parameterized by & Kq with 
the same slope created as 8 moves from 0° to ±90°. As derived from 
equations 4 and 5 above, the following formula is used to create the plane 
plot of Figure 6: 

20 (cp 23 /3Tc) +2jq 3 -/3« (<p 12 /2tO +* 12 

The ambiguity resolution algorithm is as follows. The phase 
difference between the antennas is measured as (f\ 2 and (p^. On the 
interferometer ambiguity plane of Figure 6, find the line represented by K l2 

25 and (ambiguity factors) that has the minimal distance, which is the actual 
probability distance, to (ft 2 and (fo. After solving the ambiguity, the present 
invention computes the exact direction of arrival using the full span of the 
airay i.e. 2.5A, which gives a phase gain (PG) of 57L 

One could select ratios of N:N + 1 of A/2, where N is the 

30 number of antennas in the array. This ratio of N:N + 1 of Wl may be 

selected to be as large as possible in order to improve the phase gain. But as 
N gets larger, the distance between the ambiguity lines gets smaller and the 
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algorithm becomes vulnerable to an error in the proper selection of K, 2 and 
Krj. The result will be a gross error. 

If the measurements distribution of (p is gaussian, then in order 
to compute the gross error probability, the "tail" - of the gaussian distribution 
5 from the distance between the ambiguity lines to "infinity" must be computed 
since the phase is folding after ±m This computation described above is only 
an appro xim a ti on and is almost true for small values of cr (cp) which is true at 
high signal to noise ratio (SNR) or high signal to interference ratio. For white 
gaussian noise, a ((p) could be approximated by IN SNR where SNR is the 
10 signal to noise ratio. 

Figure 7 is the ambiguity plane for the ratios 5:6 of 2/2. The 
distance between the ambiguity lines became less than half compared to the 
ratio of 2:3 which practically increased the gross error to be impractical. 
Thus, the solution is to add more antennas to the array. Adding more 
15 antennas at different ratios spreads out the ambiguity lines in hyper space 
(multi-dimensional space) with dimension of N-l (N - number of antennas in 
the array) and increases the distance between them. In the preferred 
embodiment of the present invention, it was concluded that in an urban 
environment six antennas are required in a linear array resulting in a five 
20 dimensional plot' of Figure 6 (which is incapable of being shown in a 

drawing). The following formula are used to create the plane plot of Figure 7 
which is an example of moving the antenna spacing closer together: 



25 



SIN(Q) = (k/D 12 ) * (<p 12 +27C*i 2 ) / (2tc) = (<p 12 /57c) +2K 12 /5 



SIN(6) = (k/D 23 ) *(<P 23 +2tzK 22 ) / (2ti) = (<p 23 /67t) +K 22 /3 
(cp 23 /67r) +K 22 /3= (cp 12 /57i) +2iC 12 /5 

Direction Finding Algorithm 
30 The full direction finding algorithm relates to measurements of 

a linear antenna array having six dipole antennas per side of a four-sided 
array with a single omni reference antenna (see Figure 3). The calculations 
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are performed on a six antenna subarray after the sub-array is selected to be 
closest to the SMR remote mobile unit as described above in conjunction with 
Table 1. The full interferometer formula for a six antennas sub-array is: 



<pl 












«1 


<p 2 


_ 2tz 














^3 


SIN{Q) +27t 




+ 




<p 4 


X 






*4 




c 4 


<p 5 




D 5 








E 5 



10 where 

CR is the electrical phase measured between the antennas, 
Dj is the distances between the antennas 
X is the wavelength 
0 is the direction of arrival (DOA) 
15 ' Kj is the ambiguity factor, integer value such that, 

-n < 2n Di * SIN (6) / X + 2k K { < n 
Ej is the additive noise (assumed to be gaussian). . 
To solve the equation is (6) above, linear array calculations are 
performed to solve the following values: . 



20 
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E 4 /27T 






e s /27i 



X=0 .5 SXiV(e) 
9 = d X+k+E — <P-k = d X+ E. 



(7) 



The results of these calculati ons are as follows: 

the vector elements (p are limited by: | -0.5< (p s < 0.5 | ; 
(p is the measurements vector and d is known; 
5 -0.5 < X < 0.5; X is unknown; 

k is unknown; its elements are integer, limited by the elements 

- of d 

E is the noise vector that is unknown. 
In the calculations, it is assumed that the noise is additive such as: 

10 



15 
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(8) 



1° -i=J 



tffej}* 0; i 6(1,2,3,4, 5) 



2 CT 2 i=j 
i-j-1 
i-j+l 
0 otherwise 



-a 2 
-a 2 
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R-CO 



2-10 0 0 
-12 -10 0 
0-12-10 
0 0-12-1 
0 0 0-12 
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Where : 

<j>, is the electrical phase measured at antenna i. 
. <R = <k - <ta 

25 The solution to this equation is to search among all K e N, the 

specific K that minimizes a weighting function in a weighted least square 
procedure (WLS). For every KeNwe compute: 



30 
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The error will be: 



5 



0^) * / W/fTp-l \, 

10 Weighting function: 

( H ) C(£) =£ T (£)i?'- 1 £(£) - 



15 



20 

Where A is 5 X 5 matrix. The rank of A is only 4. In order to 
' reduce the amount of real time computation, the following computations are 
performed and the results stored . 

25 (12) A=F T F 

A =A v A 0 Aj 

30 Since the rank of A is 4, one of the eigenvalues and thus the 

. corresponding eigenvector are 0. Matrices 4X5 and 5X4 are created and 
designated by a — over the top of the matrix identifier and thus: 
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(13) 



A = F Tp 



Using these values in the weighting function, produces: 

10 (14) , ,r. / , 

C(£) = (2-4 F*F (<£-£) 

<P=F<£ 

K-FK 
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It 











The full algorithm thus becomes: 
Compute in off-line: 



20 



1 . Compute F 

2 . For all, KEN compute n= 

3 . • Compute the matrix ff. 



25 



Compute in real time: 
1. Compute: 



<P=F<p 



30 



2. For-all n e N, compute: . 
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C{K) = cp-£ 



3. Find the best three (3) K that mostly minimize C(K). This procedure is 
5 done for each hop. 

4. Compute X(K) for those K selected 

5. Compute a histogram from all the X(K), using for each DOA its 
distance from the lines (weight) in inverse as weight 

6. Filter the histogram by passing a boxcar low-pass filter of 3° over the 
10 histogram in order to cancel small peaks in the neighborhood of a real 

peak 

7. Find the highest peak and the second to the highest peak. 

8. Report the results. ■ ■- . - - 

15 Trianaula fion Process 

The input to the location algorithm are reported measured 
directions from several base stations. Each base station reports two directions. 
. Each direction is reported with its peak heights which is an indication of its 
validity. The location algorithm projects the direction of the incoming signal 
20 on a gnomonic plane tangent to earth, it filters out the false readings, and it 
computes the most probable location of the SMR transmitter. 

■* 

Cassini UTM Projection . 
The solution of triangulation on the earth ellipsoid is a highly 

25 non-linear problem. In order to simplify the problem, we transfer the 

observed data to a plane called the gnomonic plane as shown in Figure 9. 
The geographic location of the base station on earth is transferred to an X-Y 
position on the gnomonic plane. Each measured bearing is also transferred to 
the gnomonic plane with a small adjustment. The estimation process is done 

30 on the plane. The estimated location of the SMR remote mobile unit is not 
transferred back to an earth ellipsoid since the same plane is used for graphic 
display of maps using a Microsoft Windows program or the like. 
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The Cassini UTM projection is widely used for mapping. Its 
advantages are that the mapping are in the Cartesian coordinate system such 
that the E-W and N-S lines are orthogonal. The projection maintains 
approximately both the directions related to the north and the distances at 
5 radius less than 50 M The projection is very convenient for graphic display 
such as in the Microsoft Windows format. 

The main distortion of the UTM projection is caused by the 
convergence of the meridian toward the north. There are two corrections that 
have to be made: 

10 1. Correction of the map north to the geographical north, which is 

a bias correction; and 
2. Small correction, of the DOA (direction of arrival) of the signal 
according to the location of the base station. 
In actual observations, the correction was actually measured 
15 together with the calibration and measurement of the pointing of the direction 
finding array toward the north. The calibration was performed relative to 
known reference point locations in the covered area. To improve system 
accuracy, a reference transmitter having a known exact location may be used 
to correct the system parameters. 

20 

Direction Finding Fixing on a Plane 
The fixing process input are the two DOA's received from all 
the base stations. Those skilled in the art will readily recognize that the 
present invention can operate with more or less DOA's from the base stations. 
25 There are several tests performed on the results received from the base 
stations (BS): 

If the peaks are too low, both measurements are ignored 

If the ratio between peaks is larger than 4, the second peak is 

ignored 

30 The CS (Central System) computes the fix by using all the combinations from 
the DOA received from the base stations. Note that within one message there 
could be several DO As. 



